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Abstract 

The presence of a sizeable CP-violating phase in Bg-Bg mixing would be an unambiguous signal 
of physics beyond the Standard Model. We analyse various possibilities to detect such a new phase 
considering both tagged and untagged decays. The effects of a sizeable width difference AF between 
the Bg mass eigenstates, on which the untagged analyses rely, are included in all formulae. A novel 
method to find this phase from simple measurements of lifetimes and branching ratios in untagged 
decays is proposed. This method does not involve two-exponential fits, which require much larger 
statistics. For the tagged decays, an outstanding role is played by the observables of the time- 
dependent angular distribution of the Bg J /+/^] K'^K^] decay products. We list the 
formulae needed for the angular analysis in the presence of both a new CP-violating phase and a 
sizeable AF, and propose methods to remove a remaining discrete ambiguity in the new phase. This 
phase can therefore be determined in an unambiguous way. 
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1 Introduction 

The rich phenomenology of non-leptonic B decays offers various strategies to explore the phase 
structure of the Cabibbo-Kobayashi-Maskawa (CKM) matrix [|T]] and to search for manifestations 
of physics beyond the Standard Model [0]. Concerning the latter aspect, CP violation in Bg-Bs 
mixing is a prime candidate for the discovery of non-standard physics. In the first place the Bg-Bg 
mixing amplitude is a highly CKM- suppressed loop-induced fourth order weak interaction process 
and therefore very sensitive to new physics. Moreover in the Standard Model the mixing-induced 
CP asymmetries in the dominant Bg decay modes practically vanish, because they are governed by 
the tiny phase axg{—VthV^g/ {VchV*g)). It does not take much new physics to change this prediction: 
already a fourth fermion generation^ can easily lead to a sizeable new CP-violating phase in Bg- 
Bg mixing It is further possible that there are new flavour-changing interactions which do not 
stem from the Higgs-Yukawa sector. The phases of these couplings are not related to the phases 
of the CKM elements and therefore induce extra CP violation. An example is provided by generic 
supersymmetric models in which new flavour-changing couplings come from off-diagonal elements 
of the squark mass matrix [|]]. While such new contributions are likely to affect also B^-Bd mixing, 
they appear in the Bd system as a correction to a non-zero Standard Model prediction for the mixing- 
induced CP asymmetry, which involves the poorly known phase (3 = a,rg(—VtbV^*d/ {VcbV*d)). To 
extract the new physics here additional information on the unitarity triangle must be used. In the Bg 
system, however, the new physics contribution is a correction to essentially zero [^. 

Indeed, the discovery of new physics through a non-standard CP- violating phase in Bg-Bg mix- 
ing may be achievable before the LHCb/BTeV era, in Run-II of the Fermilab Tevatron. 

-meson decays into final CP eigenstates that are caused by 6 — » ccs quark-level transitions 
such as Bg — > DfDj, J/ipi]^'^ or J/ilxp, are especially interesting [0-0]. The r] and rj' mesons in 
Bg J/il)ri^''> can be detected through 77 — > 77 and rj' p°7, tx^tx^t], or through 77 — > tt+tt^tt^ 
[|T0|]. These modes require photon detection. In the case of i?s J/il)[—^ /+/^] </)[^ fC+J^^], which 
is particularly interesting for 5-physics experiments at hadron machines because of its nice experi- 
mental signature, the final state is an admixture of different CP eigenstates. In order to disentangle 
them, an angular analysis has to be performed [[Tl], |12|]. Experimental attention is also devoted to 
three-body final states [|T3]]. i?s -meson decays triggered by the quark decay h cuJcan likewise 
access a CP-specific final state, e.g. via Bg -D^p^ [-^ K^K^]Ks, with a likewise negligibly small 
CP- violating phase in the Standard Model. The key point here is that there are many different decay 
modes which all contain the same information on the pursued new CP-violating phase 0. Further- 
more, additional information on (f) can be gained from analyses that require no tagging. Untagged 
studies determine | cos0| and are superior to tagged analyses in terms of efficiency, acceptance and 
purity. However, they require a sizeable width difference | Ar| between the Bg mass eigenstates. On 
the other hand, from tagged analyses (such as CP asymmetries) sin (p can be extracted, if the rapid 
Bg-Bg oscillation can be resolved. Both avenues should be pursued and their results combined. 



' This scenario is still possible, though somewhat disfavoured by electroweak precision data [^. 
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because they measure the same fundamental quantities. 

If we denote the Standard Model and the new physics contributions to the Bg-Bg mixing am- 
plitude with and Snp, respectively, then the measurement of the mass difference Am in the Bg 
system determines I^sm + •S'npI- The knowledge of both Am and the Bg-Bg mixing phase (j) then 
allows to solve for both the magnitude and phase of Snp- Information on is especially valuable, if 
I^smI and |S'np| are comparable in size and Am agrees within a factor of 2 or 3 with the Standard 
Model prediction. 

The purpose of this paper is twofold: we first identify useful measurements and show how the 
information from different decay modes and different observables can be combined in pursuit of a 
statistically significant "smoking gun" of new physics. Second we show how the Bg-Bg mixing 
phase can be identified unambiguously, without discrete ambiguities. The outline is as follows: after 
setting up our notation in Section |^ we consider untagged Bg decays and discuss various methods 
to determine | cos0| in Section |[ Tagged Bg decays are discussed in Section ^, whereas Section |] 
shows how to resolve the discrete ambiguity in 0. Finally, we conclude in Section |^. 



2 Preliminaries 

In this section we define the various quantities entering the time evolution of Bg mesons and their 
decay amplitudes. We closely follow the notation of the B4B4R-Book [|T]]. Some of the discussed 
quantities depend on phase conventions and enter physical observables in phase-independent combi- 
nations [[T4|]. Since this feature is well understood and extensively discussed in the standard review 
articles [|T]], we here fix some of these phases for convenience and only briefly touch this issue where 
necessary. 

We choose the following convention for the CP transformation of meson states and quark cur- 
rents:! 

CP\ Bg) = -\Bg), CPq^i^bL (CP)-' = -hYQL. (1) 
Hence the CP eigenstates are 

\Br") = ^{\Bg)-\Bg)), and | fi^ ) = {\ Bg) + \Bg)) . (2) 
The time evolution of the Bg-Bg system is governed by a Schrodinger equation: 

4J\m)-(M-^l)(\m) (3) 



dt\\Bg{t))J V 2j\\Bg{t)) 

with the mass matrix M = and the decay matrix F = Here | Bg{t) ) denotes the state 
of a meson produced as a Bg at time t = 0, with an analogous definition for | Bg{t) ). The off- 
diagonal elements M12 = M21 and F12 = F21 correspond to Bg-Bg mixing. In the Standard Model 
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s u,c,t b 

Figure 1 : Bg-Bg mixing in the Standard Model. 



the leading contributions to M12 and stem from the box diagram in Fig. |I]; T12 originates 
from the real final states into which both Bg and Bg can decay. It receives contributions from box 
diagrams with light u and c quarks. Since is dominated by CKM-favoured tree-level decays, 
it is practically insensitive to new physics. On the other hand, M12 is almost completely induced 
by short-distance physics. Within the Standard Model the top quarks in Fig. |T] give the dominant 
contribution to Bg-Bg mixing. This contribution is suppressed by four powers of the weak coupling 
constant and two powers of \Vtg\ ~ 0.04. Hence new physics can easily compete with the Standard 
Model and possibly even dominate Mu- If the non-standard contributions to Mu are unrelated to 
the CKM mechanism of the three-generation Standard Model, they will affect the mixing phase 



Hi 



argMi2. 



With our convention ([I|) the Standard Model prediction is 0m = arg(Vtf,V^ 



*\2 
ts) ■ 



The mass eigenstates at time t = 0, \ Bl) and | Bh ), are linear combinations of\Bg) and \ Bg): 
lighter eigenstate: \ Bl) = p\ Bg) + q\ Bg) 



heavier eigenstate: | Bh ) = p\ Bg) — q\ Bg), with |pp + |gp = 1 

We denote the masses and widths of the two eigenstates with Ml,h and Tl^h and define 

TB. 2 ' 



Am 



M, 



H 



Ml, 



Ar 



(4) 



(5) 



While Am > by definition, AF can have either sign. Our sign convention is such that AF > in 
the Standard Model. By examining the eigenvalue problem of M — iF/2 we find that the experimental 
information Am ^ F model-independently implies IF12I ^ |Mi2|. By expanding the eigenvalues 
and q/p in F12/M12, we find 



Am = 2|M- 



12 



AF = 2 IF12I cos I 



and 



1 
P 



—e 



1 - 



(6) 



^metric = (1, -1, -1,-1) 
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Here the phase is defined as 



Ml 



12 



12 



M- 



12 



12 



3«0 



In @ we have kept a correction in the small parameter 

ri2 



M- 



12 



sm I 



(7) 



(8) 



but neglected all terms of order r^2/^i2 do so throughout this paper. Since a can hardly exceed 
0.01 we will likewise set it to zero in our studies of Bg decays into CP eigenstates and only briefly 
discuss a non-zero a in sect. p4[ 

The phase is physical and convention-independent; if = 0, CP violation in mixing vanishes. 
In the Standard Model = (pM — arg(— ri2) is tiny, of order 1%. This is caused by two effects: 
first, is dominated by the decay b ccs and (KbV^*)^ is close to the Bg-Bg mixing phase 
aigCVtbV^iy . Second, the small correction to arg(— ri2) involving VubV*g is further suppressed by a 
factor of mil ml. In the search for a sizeable new physics contribution to these doubly Cabibbo- 
suppressed terms proportional to VuhV*g can safely be neglected, as we do throughout this paper. 

For a Bg decay into some final state /, we introduce the | Ai?| = 1 matrix elements 



As = {f\Bg) 
The key quantity for CP violation reads 



A 



/ 



and 



pAf- 



Af = {f\Bg). 



(9) 



The time evolution formulae and the expressions for the CP asymmetries in the forthcoming sections 
can be conveniently expressed in terms of 

^'r^^ = ^^^, ^c^? = -:^^ and ^Ar = -:^^. (10) 



i + \Xf\ i + IA/i 1 + 

If / is a CP eigenstate, CP| / ) = ±| / ), then A^^. 7^ or A"^^ 7^ signals CP violation: a non- 
vanishing AqI, implies \Af \ 7^ \Af\, meaning direct CP violation; A!^^ measures mixing-induced 
CP violation in the interference of Bg ^ f and Bg f. The third quantity, ^Ar^ plays a role, if 
Ar is sizeable. The three quantities obey the relation 



A 



dir 
CP 



+ 



/) mix 



+ 1^, 



AF 



1. 



The time-dependent decay rate T{Bg{t) f) of an initially tagged Bg into some final state / is 
defined as 



1 dN{Bg{t)^f) 



B 



dt 



(11) 
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Quark decay 



CCS 



cud 



Hadronic decay 



Bs 
Bs 
Bs 
Bs 
Bs 
Bs 
Bs 
Bs 
B. 



n(*)+ n(*)- 

s s 

/}(*)+£)(*)- or £)(*)o_d(*)o 



Bs KsD^*^^ (pKs, p^Ks, KK or tt+tt-] 



Remarks 



CP-odd final state 
CP-odd final state 
D^D- is CP-even 
non-spectator decays, 
DD is CP-even 



Table 1: Some CKM-favoured Bs decay modes into CP-specific final states. Here, i}) represents J/ip 
or il){2S). Decays into two vector particles or into three-body final states with one or more vector 
particles require an angular analysis to separate the CP-even from the CP-odd component. The final 
states DfD*^ are dominantly CP-even [ p^ (see sect. |3|). 



Here Bs{t) represents a meson at proper time t tagged as a 5^ at t = 0; dN{Bs{t) /) denotes the 
number of decays of Bs{t) into the final state / occurring within the time interval [t, t + dt]; Nb is 
the total number of Bs's produced at time t = 0. An analogous definition holds for T(Bs{t) — * /). 
By solving the Schrodinger equation (U) using @, we can find these decay rates [fT^]: 

,2l + |A/|' 



TiBsit) -^f)=Aff \A 



f\ 



-rt 



Art Art 

cosh h ^cp cos(Amt) + Aaf sinh h ^cp sin (Ami) 



,(12) 



T{Bs{t)-^f)=^f\Af\ 



|A/r 



'1 + a) e 



-rt 



X 



cosh 



Art 



^cp cos(Amt) + sinh 



Art 



A^l^ sin(Amt) 



(13) 



Here A/} is a time-independent normalization factor. 

A promising testing ground for new physics contributions to (pM are decays into CP eigenstates 
triggered by the quark decay b ccs. Table [1] summarizes such CP-specific Bs decay modes. To 
estimate the size of the small Standard Model predictions consider first the decay amplitudes [[TtI]: 

A2 \ 



Af.Af oc 



1 + 



(14) 
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Hence the weak phase factor e^"', which is associated with the quantity OpC*^, is strongly Cabibbo 



suppressed by two powers of the Wolfenstein parameter A ~ \Vus\ — 0.22 [p^. The "penguin 
parameter" OpC*^ measures - sloppily speaking - the ratio of penguin- to tree-diagram-like topologies 
and is loop-suppressed. Since new-physics contributions to these decay amplitudes have to compete 
with a tree diagram, they are not expected to play a significant role. A detailed discussion for a 
left-right-symmetric model can be found in [Q]. Since we are interested in large "smoking gun" 
new physics effects in Bg-Bg mixing, we account for the Standard Model contributions within 
the leading order of A and set \Af\ = \Af\, neglecting direct CP violation. With the weak phase 
(pecs = arg(KfeKs) oils then finds 

^ = -7]fe^'^'^-^^. (15) 

Here r]f denotes the CP parity of/: CP\ f) = r]f \ f). In Table [l| we also included decay modes 
driven by the quark level decay b cud. The weak phase of these modes involves the phases of the 
K and D decay amplitudes into CP eigenstates. The phases combine to argCVcbV*^) + aTg{Vudy*g) + 
aLrg{VusV*g) = aTg{VcbV*g), i.e. the same result as for b ccs. With @ and ((T3|) A/ reads 

Xf = -^ = Vfe-'t (16) 

Here we have identified the phase arg{r]f\f) = (pM — with the phase </> defined in ([?[). This is 
possible, because arg(— ri2) = 20ccs + 0(A^) and we neglect the Cabibbo-suppressed contributions. 
The Standard Model contribution to = 0sm + 0np equals 0sm = —2r]\'^. Here t] is the Wolfenstein 
parameter measuring the height of the unitarity triangle. Since our focus is a sizeable new physics 
contribution 0np, we can safely neglect 0sm and identify with 0np in the following. That is, we 
neglect terms of order A^ and higher. Using ( [T^ ) the quantities in (|TD|) simplify to 

Acp = 0, ^cp='7/sin0 and Aaf = -Vf coscp. (17) 

The corrections to (|17]) from penguin effects can be found in [jlTp. We next specify to the PDG 
phase convention for the CKM matrix [|T^], in which arg(V^feV^*) = 0{\^). Then we can set (pczs to 
zero and identify 

<pM = 0- 

With this convention the mass eigenstates can be expressed as 

I bl ) = 1 5r ) — Y~ ' ^ + 

1 _ i0 1 _|_ 10 

\Bu) = - ^— I 5r ) + I Bf') + 0{a) . (18) 

Whenever we use Bf^^ and 5°'*'^ we implicitly refer to this phase convention. If formulae involving 
j^even ^odd used to coustraiu models with an extended quark sector, the phase convention 
used for the enlarged CKM matrix must likewise be chosen such that arg(VcfeV^* ) ~ 0. 
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3 Untagged Studies 
3.1 Time Evolution 

Whereas the width difference AF is negligibly small in the 3^ system, it can be sizeable for Bg 
mesons. This has the consequence that the untagged Bg data sample bears information on CP vio- 
lation [^. Further the width difference itself is sensitive to the Bg-Bg mixing phase (j) [|2l|], as we 
can see from @. 

When i?s's and .B^'s are produced in equal numbers, the untagged decay rate for the decay B™ 
f reads 

T[f,t] = T{Bg{t)^f) + T(Bg{t)^f) 



Aff 



~rrt 



\{f\BL)\' + e-'"'\{f\BH)\' +Oia). 



(19) 



1 + IA, 



ATt ATt 
cosh — — + sinh — — ^Ar \ + 0{a). (20) 



Here the second expression is simply obtained by adding ((T^) and ([T3|). In ( |T9| ) the same result is 
expressed in terms of the mass eigenstates and nicely exhibits how the decay is governed by two 



exponentials. Using dTTj ) we can relate the overall normalization to the branching ratio: 

1 r'^ 

2 Jo 



Br[f] 



dtT[f,t] 



f 



A, 



1 + |A/ 



F + ^Ar Ar/2 
F2 - (AF/2)2 



+ 0{a). 



(21) 



(22) 



Conforming with [jT^ we have normalized the event counting to Nb + N^ = 2Nb, so that Br[aU] 
1. Using (P3) we rewrite (EO) as 



T[f,t]=2Br[f] 



F + ^Ar AF/2 



, ATt . , AFt , ■ 

cosh h smh ^Ar 



+ 0{a) 



(23) 



Now ( p3| ) is our master equation for the time evolution of the decay of an untagged Bg sample. If 
F = 1/ tb^ is known, one could perform a two-parameter fit of the decay distribution to (121) and 



determine AF and ^Ar- The latter determines cj) through ([T7[), if / is a CP eigenstate from a CKM- 
favoured decay. In practice, however, most data come from short times with AF t <^ 1, and one is 
only sensitive to the product AF • ^Ar^ 



T[f,t] = 2Br[f]Te 



-rt 



\ Ar , 

1 + — ^Ar 



t 



+ 0[{ATty 



(24) 



We return to this point in sect. 



3.3, 
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3.2 The Width Difference AF and Branching Ratios 



The mass matrix M12 and the decay matrix provide three rephasing invariant quantities: IM12I, 
|ri2| and the relative phase 0. In @ we have related the two observables Am and AF to |Mi2|, 
|ri2| and (p. Interestingly, it is possible to find a third observable, which determines |ri2| and thus 
encodes additional information. We define 



AFcp = 2|r 



12 



2 E inBs 



/cP+) - TiB, ^ /, 



CP- 



(25) 



Here X^z represents the final states containing a (c, c) pair, which constitute the dominant contribu- 
tion to AFcp stemming from the decay b ccs. In ( P5[ ) we have decomposed any final state / into 

/cp+ ) + I /cp- ) and defined 



its CP-even and CP-odd component, | / ) 

T{B,-. fcp±) =Aff\{fcp± 



B. 



l(/< 



CP± 



\Bg 



\if\Bs 



A// is the usual normalization factor originating from the phase-space integration. In order to prove 
the second equality in (P5|) we start from the definition of ri2: 



T^, = Y.^ff{Bs\f){f\B.) = -Y.^ff[{Bs\f){f\Bs) + {Bs\f){f\B,) 
f ^ / 



(26) 



In the second equation we have paired the final state | / ) with its CP conjugate | / ) = —CP\ / ) . In 
the next step we trade / for /cp+ and /cp_ and use the CP transformation 



where 



{fcp±\Bs) = Te'"'^--{fcp±\B,), 
a,Tg{VcbV*g) is the phase of the b ccs decay amplitude, which dominates ri2. Then 



1) becomes 



12 



E ^ff [i(/cp+is.)r-i(/cp-is. 



E inBs 

fax-. 



/cP+) - T{B, ^ /cP-)] 



(27) 



Interference terms involving both ( /cp+ \Bs) and ( /cp- \Bs) drop out when summing the two terms 
{Bs\f){ f \Bs) and ( Bg \f) { f \Bs). In ( P7[ ) both sides of the equation are rephasing-invariant. An 
explicit calculation of ri2 reveals that the overall sign of the LHS of ( p7| ) is positive, which completes 
the proof of (^. 

Loosely speaking, AFcp is measured by counting the CP-even and CP-odd double-charm fi- 
nal states in Bg decays. We specify this statement in the following and relate AFcp to mea- 
sured observables in sect. p.3.2| . Our formulae become more transparent if we adopt the stan- 
dard phase convention with aigiVcbV*^) ~ and use the CP-eigenstates defined in (@). With 
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\Bs) = {\ 5f ™ ) + I Bf'^))/V2 one easily finds from (^: 

AFcp = 2|ri2| = r (Ef™) - r(fif (28) 

Here the RHS refers to the total widths of the CP-even and CP-odd Bg eigenstates. We stress 
that the possibility to relate |ri2| to a measurable quantity in ( ^51 ) crucially depends on the fact 
that is dominated by a single weak phase. For instance, the final state K^K^ is triggered 
by 6 ^ uus and involves a weak phase different from h ccs. Although K^K' is CP-even, 
the decay B^^'^ K^K~ is possible. An inclusion of such CKM-suppressed modes into 



would add interference terms that spoil the relation to measured quantities. The omission of these 



contributions to induces a theoretical uncertainty of order 5% on (P8|). 

In the Standard Model the mass eigenstates in (jTSp coincide with the CP eigenstates (with B^ = 
^even-j AFsm = AFcp. The effect of a non-zero B^-Bg mixing phase reduces AF: 

AF = AFcpcos0, (29) 

while AFcp = 2|Fi2| is not sensitive to new physics. From the calculated F12 we can predict to 
which extent F {Bf^") exceeds F (^-8°'*'^^ and this result does not change with the presence of a 
non-zero 0. 

The theoretical prediction for AFcp is known to next-to-leading order in both and 
the QCD coupling [ ^3| ] . It reads 

\ 2 



AFcp / I 



Bs 



245 MeV 



(0.234 ± 0.035) Bs - 0.080 ± 0.020 ] . (30) 



Here the coefficient of Bs has been updated to mb{mb) +ms(mb) = 4.3 GeV (in the MS scheme) and 
is the Bg meson decay constant. Recently the KEK-Hiroshima group succeeded in calculating 
in an unquenched lattice QCD calculation with two dynamical fermions [E?p. The result is /b^ = 



(245 ± 30) MeV. Bs parametrizes the relevant hadronic matrix element, with i?s = 1 in the vacuum 
saturation approximation. A recent quenched lattice calculation has yielded Bs = 0.87 ± 0.09 [p 



for the MS scheme. A similar result has been found in [p6|]. This analysis, however, calculates 
AF after normalizing ( PPj ) to the measured mass difference in the Bd-Bd system. This method 
involves \ Vtd\, which is obtained from a global CKM fit and thereby relies on the Standard Model. 
Since the target of our analysis is new physics, we cannot use the numerical prediction for AF 
of [E^. At present, studies oi Bs are a new topic in lattice calculations and we can expect substantial 



improvements within the next few years. With these numbers one finds from (pQ): 

AFcp 

= 0.12 ±0.06. (31) 

Here we have conservatively added the errors from the two lattice quantities linearly. 

Since AFcp is unaffected by new physics and AFcp > 0, several facts hold beyond the Standard 
Model: i) There are more CP-even than CP-odd final states in Bg decays, ii) The shorter-lived mass 
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eigenstate is always the one with the larger CP-even component in (|rS|). Its branching ratio into a 
CP-even final state /cp+ exceeds the branching ratio of the longer- lived mass eigenstate into /cp+, if 
the weak phase of the decay amplitude is close to arg KfeKs- "^os > has a shorter lifetime 
than Bh, while for coscp < the situation is the opposite pl|]. iii) Measurements based on the 
comparison of branching ratios into CP-specific final states determine AFcp rather than AF. Such 



an analysis has recently been performed by the ALEPH collaboration [27]. ALEPH has measured 



2fir[DW+L)W"] = 0.26 



+0.30 
-0.15 



(32) 



and related it to AFcp. For this the following theoretical input has been used [|T6|]: 



i) In the heavy quark limit rric ^ oo and neglecting certain terms of order l/N^ (where Nc = 3 
is the number of colours) the decay B°^'^ DfD*"^ is forbidden. Hence in this limit the final 
state in 5™ DfD*^ is CP-even. Further in 5™ D*+D*- the final state is in an S-wave. 



ii) In the Shifman-Voloshin (SV) limit mc — oo with — 2mc PSQ, AFcp is saturated 
by F(fi™ £)(*)+£)(*)-). With i) this implies that in the considered limit the width of 5f ^ 
vanishes. For Nc ^ oo and in the SV limit, 2F(i?™ D'f^^Dl*^^) further equals the parton 



model result for AFcp (quark-hadron duality). 
Identifying F(fi|™" £)(*)+£)(*)-) ^ AFcp and T{Bf'^ 



2 5r[DW+L)W-] ^ AF 



CP 



1 



COS( 



2F, 



+ 



1 



cos ( 



2F 



H 



nM+ni*)- 

s s 

AFcp 



we find: 



(33) 



Thus the measurement in ( |32D is compatible with the theoretical prediction in (pT]). For = 0, 
the expression used in Ref. [^7p, in which the Standard Model scenario has been considered, is 
recovered. The term in square brackets accounts for the fact that in general the CP-even eigenstate 
I B^^'^^ ) is a superposition of | B^ ) and | Bh )■ It is straightforward to obtain (|33|): inserting ( [T8| ) into 
(jl^) expresses F[/, t] in terms of F(5f ™ /) and V{B'f'^ /). After integrating over time the 
coefficient of F(i?^™" f) is just the term in square brackets in i 



When using ( P3| ) one should be aware that the corrections to the limits i) and ii) adopted in [|T^ 
can be numerically sizeable. For instance, in the SV limit there are no multibody final states like 
D'f^DXg, which can modify (^3]). As serious would be the presence of a sizeable CP-odd component 
of the D'f^'^D^*^ ^ final state, since it would be added with the wrong sign to AFcp in (|33|). A method 
to control the corrections to the SV limit experimentally is proposed in sect. p.3.2| . We further verify 
from ( p3l ) that the measurement of Br[D'f^'^D^/^^] determines AFcp. Its sensitivity to the new 
physics phase is suppressed by another factor of AF /F and is irrelevant in view of the theoretical 
uncertainties. 



3.3 Determination of AF and | cos 0| 

There are two generic ways to obtain information on AF: 
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i) The measurement of the Bs lifetime in two decay modes 5™ /i and 5™ f2 with 
^Ar(/i)7^^Ar(/2). 

ii) The fit of the decay distribution of _B™ — » / to the two-exponential formula in (|23|). 



As first observed in [PT|], the two methods are differently affected by a new physics phase (p ^ 0. 
Thus by combining the results of methods i) and ii) one can determine 0. In this section we consider 
two classes of decays: 

• flavour- specific decays, which are characterized by Af = implying = 0. Examples are 

Bs D-TI+ and 5, Xi+ue, 

• the CP-specific decays of Table |I], with = —Vf cos (p. 



In both cases the time evolution of the untagged sample in ( p3| ) is not sensitive to the sign of AF (or, 
equivalently, of cos 0). For the CP-specific decays of Table [I| this can be seen by noticing that 

• u ^rt , ^. ^ \AT\t 

Aat smn ^ = ~Vf \ cos0| smn — - — . 

Here we have used the fact that AF and cos always have the same sign, because AFcp > 0. Hence 
the untagged studies discussed here in sect. ^3] can only determine | cos0| and therefore lead to a 
four-fold ambiguity in 0. The sign ambiguity in cos reflects the fact that from the untagged time 
evolution in (531) one cannot distinguish, whether the heavier or the lighter eigenstate has the shorter 



lifetime (however, see sect. 5). 

In order to experimentally establish a non-zero AF from the time evolution in (P3[) one needs 



sufficient statistics to resolve the deviation from a single-exponential decay law, see (P4|). As long as 
we are only sensitive to terms linear in AF t and AF/F, we can only determine Aaf AF from (P^. 
.Aaf Ar vanishes for flavour- specific decays and equals —7]fAT cos0 for CP-specific final states. 
Hence from the time evolution alone one can only determine AF cos0 in the first experimental 
stage. This determination is discussed in sect. p.3.1| . Once the statistical accuracy is high enough to 
resolve terms of order (AF)^, one can determine both |AF| and | cos0|. Fortunately, the additional 
information from branching ratios can be used to find |AF| and | cos0| without resolving quadratic 



terms in AF. The determination of lAFl and I cos 01 is discussed in sect. 3.3.2 



3.3.1 Determination of F and AF cos 



Lifetimes are conventionally measured by fitting the decay distribution to a single exponential. Con- 
sider a decay which is governed by two exponentials. 



r[f,t] + T[f,t] 



Ae 



-Tr.t 



+ Be 



-rt 



ATt ATt 
{A + B) cosh V{B - A) sinh 



(34) 
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but fitted to a single exponential 



rft 



F[f,t] = Tfe~'f 



(35) 



In ( P^ we have averaged over / and its CP-conjugate /. Of course the coefficients depend on the 
final state: A = A(f), B = B{f). A maximum likelihood fit of ( p5| ) converges to [ pT| ] 



A/Tl + B/Th 

Ain+B/vi- 



We expand this to second order in AF: 



A-B^T 2AB {ATf 



r + — - — 



A + B 2 {A + By r 



(36) 



+ o 



'{ATf 
Y2 



(37) 



In flavour-specific decays we have A = B (see (P3[)). We see from ( P7| ) that here a single-exponential 
fit determines T up to corrections of order AF^/F^. 

Alternatively, one can use further theoretical input and exploit that F^^/F^^ = 1 + from 



heavy quark symmetry [[22|, [29|, pOQ . This relation can therefore be used to pinpoint F in terms of 
the well-measured B^ lifetime. New physics in the standard penguin coefficients of the effective 
AB = 1 hamiltonian only mildly affects TbJTb^ [Q- The full impact of new physics on TbJTb^, 
however, has not been studied yet. 

With ( P3| ) and (Q) we can read off A and B for the CP-specific decays of Table |T] and find 
A(/cp+)/5(/cp+) = (1 + cos0)/(l - COS0) and A(/cp-)/5(/cp-) = (1 - cos0)/(l + cos</)) 
for CP-even and CP-odd final states, respectively. Our key quantity for the discussion of CP-specific 
decays 5™ — > /cp is 



af: 



CP 



-r^j^ApAF = AF cosi 



AF 



CP cos 



(38) 



With this definition ( p7| ) reads for the decay rate Fcp,r;f measured in 5™ fcp' 



^CP,m = ^ + Vf 



AF' 



CP 



sm 



(ATf 
2F 



/(ATf 
1^ F2 



That is, to first order in AF, comparing the B™ lifetimes measured in a flavour-specific and a CP- 
specific final state determines AF'^p. Our result agrees with the one in [pT]|, which has found ( ^8] ) 
by expanding the time evolution in ( |^ ) and (^5]) for small AFt. Including terms of order (AF)^, 
lifetime measurements in a flavour-specific decay 5™ determine []3T|]: 



F - 



(ATT , ( W 
2F V ^2 ^ 
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This implies Ffs < F. Despite the heavy quark symmetry prediction T Bs/'^b^ — 1, a large AF leads 
to an excess of the Bg lifetime measured in 5™ — > /fs over the Bd lifetime [|T]]. From (^7]) one finds 

rcP,„ - r,. . ^ . ^) . o ((^) . (39, 

Hence for a CP-even (CP-odd) final state the quadratic corrections enlarge (diminish) the difference 
between the two measured widths. A measurement of AF'(^p at Run-II of the Tevatron seems to be 
feasible. The lifetime measurement in the decay mode B™ J li^<\> has been studied in simulations 
[l32|,|33|]. This decay mode requires an angular analysis to separate the CP-odd (P-wave) from the 



CP-even (S-wave and D-wave) components. The angular analysis is discussed in sect. |4^. With 
2 fb^^ integrated luminosity CDF expects 4000 reconstructed 5™ — > J [H^<\> events and 
a measurement of AF^p/F with an absolute error of 0.052. This simulation assumes that F — 
(AF)^/(2F) (see (^7])) will be measured from flavour-specific decays with an accuracy of 1% 
and uses the input AF^p/F = 0.15. When combining this with other modes in Table |I] and taking 
into account that an integrated luminosity of 10-20 fb~^ is within reach of an extended (up to 2006) 
Run-II, the study of AF'^p at CDF looks very promising. The LHC experiments ATLAS, CMS and 
LHCb expect to measure AF'^p/F with absolute errors between 0.012 and 0.018 for AF'^p/F = 0.15 
[|34l]. An upper bound on AF^p would be especially interesting. If the lattice calculations entering 
( |3T| ) mature and the theoretical uncertainty decreases, an upper bound on | AFqp| may show that 
7^ 0, TT through 

^ = cos^0. (40) 

Ai CP 

Note that conversely the experimental establishment of a non-zero AF'(^p immediately helps to con- 
strain models of new physics, because it excludes values of around 7r/2. This feature even holds 
true, if there is no theoretical progress in (|3T1). 

The described method to obtain AF'^p can also be used, if the sample contains a known ratio 
of CP-even and CP-odd components. This situation occurs e.g. in decays to J/V'0, if no angular 
analysis is performed or in final states, which are neither flavour-specific nor CP eigenstates. We 
discuss this case below in sect. [3.3.2| with 5™ DfD'f^^. A measurement of the Bg lifetime in 



-B™ — > J/ Tpcj) has been performed in [^, but the error is still too large to gain information on AF^p. 
Note that the comparison of the lifetimes measured in CP-even and CP-odd final states determines 
AF'f.p up to corrections of order (AF/F)^. 

3.3.2 Determination of | Ar| and | cos 4>\ 

The theoretical uncertainty in ( pT| ) dilutes the extraction of | cos0| from a measurement of AF'^p 
alone. One can bypass the theory prediction in ( 01| ) altogether by measuring both AF'^p and | AF| 
and determine | cos0| through 

AF' 

^f^ = |cos0|. (41) 
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To obtain additional information on AF and (p from the time evolution in ( P3| ) requires more statis- 
tics: the coefficient of t in (|24l), Ar^Ar/2, vanishes in flavour-specific decays and is equal to 
—r]fAT'Qp/2 in the CP-specific decays of Table [I]. Therefore the data sample must be large enough 
to be sensitive to the terms of order (AF t)^ in order to get new information on AF and </>. We now 
list three methods to determine |AF| and | cos0| separately. The theoretical uncertainty decreases 
and the required experimental statistics increases from method 1 to method 3. Hence as the collected 
data sample grows, one can work off our list downwards. The first method exploits information from 
branching ratios and needs no information from the quadratic (AF ty terms. 



Method 1: We assume that AF^p has been measured as described in sect. p.3.1| . The method 
presented now is a measurement of AFcp using the information from branching ratios. With (^0|) 
one can then find | cos0| and subsequently |AF| from In the SV limit the branching ratio 

+/}(*)-] equals AFcp/(2F) up to corrections of order AF/F, as discussed in sect. ^ 
[1T6[]. Corrections to the SV limit, however, can be sizeable. Yet we stress that one can control the 
corrections to this limit experimentally, successively arriving at a result which does not rely on the 
validity of the SV limit. For this it is of prime importance to determine the CP-odd component 
of the final states DfDl^ and D*~^D*^. We now explain how the CP-odd and CP-even component 
of any decay 5™ f corresponding to the quark level transition b — > ccs can be obtained. This 
simply requires a fit of the time evolution of the decay to a single exponential, as in (P5|). Define 
the contributions of the CP-odd and CP-even eigenstate to Bs ^ f: 



TiBf^^f) ^ Ar^|(/|i?r)p 



T{Br~^f) ^ Aff\{f\Bl 



(42) 



It is useful to define the CP-odd fraction x/ by 



T{B 



odd 



r{Bi 



{f\Bf') 



■/) \{f\B: 



even 
s 



(Jl^even) 



1-Xf 



(43) 



The time evolution (F[/, t] + T[f,t])/2 of the CP-averaged untagged decay B^ 
by a two-exponential formula: 



/, / is governed 



F[/,t]+F[/,t] ^ ^^^^ ^_r,i ^ ^^^^ ^_r,,_ 



(44) 



With (|T8|) and (jT^) one finds 



Mf) = ^\{f\BL)\' 



2 

1 + cos ( 



TiBl 



^\{7\B,)\' 
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B{f) = ^\{f\Bn\ 



1 — cos I 



ViBl 



f) + 



1 + cos I 



With we arrive at 

A{f) {l + cos<P)T{Bl 



/) + (!_ cos 0)r(s: 



odd 



V{B 



f) 



odd 



/)■ 



1 + (1 -2x/ 



COS( 



COSi 



(45) 



-. (46) 



B{f) (1 - cos 0)r(fif ™ ^ /) + (1 + COS 0)r(fi°dd ^ J) 1 _ (1 _ 2x 

In ( |45| ) and (^6[) it is crucial that we average the decay rates for i?™ ^ / and the CP-conjugate 
process fi™ J. This eliminates the interference term {Bf'^ |/)(/ l^r")' so that A{f)/B{f) 
only depends on Xf. The single exponential fit with ( P5| ) determines F/. Equations (^7|) and (| 
combine to give 



(47) 



2 (Fj - r) = (1 - 22;/) Ar COS0 = (1 - 2xf) AFcp cos^ (f) 

up to corrections of order (Ar)^/r. In order to determine xj from ( |47| ) we need AF'^p from the 
lifetime measurement in a CP-specific final state like DfDj or from the angular separation of the 
CP components in 5™ ipcp. The corrections of order (Ar)^/r to ( pT] ) can be read off from ( ^T\ ) 
with (^ as well. Expressing the result in terms of F/ and the rate Ffs measured in flavour-specific 
decays, we find 



l-2x 



f 



AF 



CP 



1-2- 



(48) 



In order to solve for F(Bf ™ ^ /) and r{Bf'^ 
Recalling (||) one finds from (g^) and (|5|): 



/) we also need the branching ratio Br[f] +Br[f] . 



Br[f]+Br[f]= T{B: 



+F(5: 



odd 



/) 



/) 



1 + COS I 
2Fr 



1 — cos ( 



+ 



COS( 



2F 



H 



1 + COS ( 



By combining 



2Tl 2Th 
1) and (B9|) we can solve for the two CP components: 



(49) 



T{B 



even 

s 



/)= F2-(AF/2)^ (Br[f]+Br[f] 



odd 

s 



f) 



2F-F/ 

l-Xf) [Br[f]+Br[f]) F+ O (AF) 
F2-(AF/2)^] {Br[f] + Br[f]) ^^.^ 



= Xf (Br[f] + Br[f]) F + O (AF) 
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From we now find the desired quantity by summing over all final states /: 

AFcP = r (Bin - r (Bf') =2\t'- (Ar/2)^1 Y: Br[f] (50) 



2r Y: Br[f]{l-2xf) 



+ o 



(51) 



It is easy to find AFcp: first determine 1 — 2xf from (^8|) for each studied decay mode, then in- 
sert the result into (jS^). The small quadratic term (Ar/2)^ = ArcpArcp/4 is negligible. This 
procedure can be performed for Br[DfD*^] and Br[D*~^D*''] to determine the corrections to the 
SV limit. In principle the CP-odd P-wave component of Br[D*'^D*^] (which vanishes in the SV 
limit) could also be obtained by an angular analysis, but this is difficult in first-generation experi- 
ments at hadron colliders, because the photon from Dg'j cannot be detected. We emphasize 
that it is not necessary to separate the D^*^^D^g*^^ final states; our method can also be applied to 
the semi-inclusive D'f^'^Dl*^^ sample, using AF'^p obtained from an angular separation of the CP 
components in i?™ — > ipcp. Further one can successively include those double-charm final states 



which vanish in the SV limit into (pOj). If we were able to reconstruct all b — > ccs final states, we 
could determine AFcp without invoking the SV limit. In practice a portion of these final states will 
be missed, but the induced error can be estimated from the corrections to the SV limit in the mea- 
sured decay modes. By comparing AFcp and AF'(-.p one finds | cos0| from (p^. The irreducible 
theoretical error of method 1 stems from the omission of CKM-suppressed decays and is of order 

2|KfeW(K6K.)| -5%. 

Method 1 is experimentally simple: at the first stage (relying on the SV limit) it amounts to 
counting the B"" decays into D^*^^ D^*^^ . A first simulation indicates that CDF will be able to 



separate the Bg decay modes into DfD^ , D*^Df and D*^D* [B^. The corrections to the SV 



limit are obtained by one-parameter fits to the time evolution of the collected double-charm data 
samples. This sample may include final states from decay modes which vanish in the SV limit, such 
as multiparticle final states. No sensitivity to (AF t)^ is needed. A further advantage is that AFcp is 
not diminished by the presence of new physics. 



Method 2: In the Standard Model the decay into a CP eigenstate /cp is governed by a single 
exponential. If a second exponential is found in the time evolution of a CKM-favoured decay fi™ — > 
/cp, this will be clear evidence of new physics [^. To this end we must resolve the time evolution 
in ( |23| ) up to order (AFt)^. At first glance this seems to require a three-parameter fit to the data, 
because F[/, t] in ( P3| ) depends on F, AF and (through ^af, see ((T^) on (f). It is possible, however, 
to choose these parameters in such a way that one of them enters F [/cp,i] at order (AF)^, with 
negligible impact. The fit parameters are F' and Y. They are chosen such that 



F[/cP+,t] = 2i?r[/cp+] FV 



-r't 



l + YT't (-1 + ^-^) +0((AF)^ 



(52) 
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Here we have considered a CP-even final state, for which a lot more data are expected than for CP- 
odd states. With ( l52| ) we have generalized the lifetime fit method described in sect. |3.3.1| to the order 
(Ar ty. A non-zero Y signals the presence of new physics. The fitted rate T' and Y are related to 
r, Ar and by 



Y 



(ATT 

4r'2 



sm 



r = r(i-r) 



cos I 



-Ar. 



(53) 



Note that for | cos0| = 1 the rate r' equals the rate of the shorter- lived mass eigenstate and the 
expansion in ( |52| ) becomes the exact single-exponential formula. After determining r' and Y we 
can solve ( |53| ) for r, Ar and </>. To this end we need the width Tf^ measured in flavour-specific 
decays. We find 



I An 



)' + rlY 



AT 



r 



r 




AT' 



CP 



2 [r' - r (1 - Y)] 



sm( 



2TVY 
\AT\ 



.(54) 



The quantity AT'^p, which we could already determine from single-exponential fits, is now found 
beyond the leading order in AT/T. By contrast, AT and | sin 0| in ( p4| ) are only determined to the 
first non- vanishing order in AT/T. 

In conclusion method 2 involves a two-parameter fit and needs sensitivity to the quadratic term 
in the time evolution. The presence of new physics can be invoked from Y and does not require 
to combine lifetime measurements in different decay modes. 



Method 3: Originally the following method has been proposed to determine | Ar| [|2y,|2T|]: The 
time evolution of a 5™ decay into a flavour-specific final state is fitted to two exponentials. This 
amounts to resolving the deviation of cosh(Ar t/2) from 1 in { ^3] ) in a two-parameter fit for r and 
I Ar I . If one adopts the same parametrization as in (^2|), r' and Y are obtained from ( |53| ) by replacing 

with it/2. The best suited flavour- specific decay modes at hadron colliders are i?™ — > D^*^'^n^, 
-B™ L)^*)^7r'F7r+7r~ and 5™ D^J*^^Xi^u. Depending on the event rate in these modes, method 
3 could be superior to method 2 in terms of statistics. On the other hand, to find the "smoking gun" 
of new physics, the |Ar| obtained must be compared to Ar'(-.p from CP-specific decays to prove 

1 cos 01 7^ 1 through (^T]). Since the two measurements are differently affected by systematic errors, 
this can be a difficult task. First upper bounds on | Ar | using method 3 have been obtained in []37|]. 

TheL3 collaboration has determined an upper bound |Ar|/r < 0.67 by fitting the time evolution 
of fully inclusive decays to two exponentials [^. This method is quadratic in Ar as well. The 
corresponding formula for the time evolution can be simply obtained from ( |3^ ) with A = T^ and 
B = Th. 
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3.4 CP Violation in Mixing and Untagged Oscillations 

In the preceding sections we have set the small parameter a in (|^) to zero. CP violation in mix- 
ing vanishes in this limit. The corresponding "wrong-sign" CP asymmetry is measured in flavour- 
specific decays and equals 

TfBJt) ^ /) - T(BJt) ^7) _ , , _ , 

Ofs = - -i ^-4^ 4 = « for Af = and L4. = Aj ■ (55) 

A special case of a^s is the semileptonic asymmetry, where / = Xi~^u. A determination of a gives 
additional information on the three physical quantities |Mi2|, |ri2| and (p characterizing Bg-Bg 
mixing. Measuring Am, AFcp, Ar'(jp and a overconstrains these quantities. 
The "right-sign" asymmetry vanishes: 

T{Bg{t) ^ /) - T(Bg{t) ^ 7) = for = and \Af\ = \Aj\. (56) 

This implies that one can measure Ofg from untagged decays. This observation was already made 
in [^9|]. It is easily verified from the sum of ( [121 ) and ( [T3t ) that to order a the time evolution of 
untagged decays exhibits oscillations governed by Am. Since a is small, one must be concerned to 
which accuracy \Af \ = \Aj\ holds in flavour- specific decays in the presence of new physics. For 
example in left-right-symmetric extensions of the Standard Model, small CP- violating corrections 
to the decay amplitude could eventually spoil this relation at the few per mille level. Further, a small 
production asymmetry e = N-^/Nb — 1 also leads to oscillations in the untagged sample. To first 
order in the small parameters a, e and — 1 one finds 

r[/,t]-r[7,t] 
r[/,t] + r[/,t] 

l^fP ~ I^tP a a + e cos(Amt) ^ — 

+ TTT TT for Af = Q and \Af\ ^ \Aj\. (57) 

\Af\^ + \Aj\^ 2 2 cosh(Art/2) ^ I /I I /i v y 



unt 
"fs 



For \Af\ = \Aj\ and e = one recovers the formula derived in p9|]. Note that the production 
asymmetry between Bg and Bg cannot completely fake the effect of a non-zero a in (^: while 
both a 7^ and e 7^ lead to oscillations, the offset from the constant term indicates new CP- 
violating physics either in Bg-Bg mixing (through a 7^ 0) or in the studied decay amplitude (through 
\Af \ 7^ l^/D- The latter effect, which is theoretically much less likely, can be tested in B^ decays 
and can therefore be disentangled from a 7^ 0. 

The ratio AFcp/F < 0.22 from (|3|) and the current experimental limit Am > 14.9 ps"^ 
imply that \a\ < 0.01. CDF expects sufficiently many reconstructed — > 0^*^"^^^ and — > 
i5^*^^7r^7r+7r^ events at Run-II after collecting 2 fb^^ of integrated luminosity to achieve a statistical 
error at the few permille level. From (|^) and (|^) we can relate a to | Ar|, Am and 0: 

I API sin0 
a = — — -. -. 

Am cosa> 
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Note, however, that the measurement of the sign of a determines the sign of sin (f). This reduces the 
four- fold ambiguity in (j) from the measurement of | cos0| to a two-fold one. It is interesting that, 
at order a, without tagging one can in principle gain information which otherwise requires tagged 
studies. Of course sin (p can be measured more directly from tagged decays, as discussed in the 
forthcoming section ^. 

4 Tagged Decays 

4.1 The CP- Violating Observables of Bs DtDj and J/^Jj rj^'^ 

For a Bs decay into a CP eigenstate / the Bg-Bg oscillations lead to the following time-dependent 
CP asymmetry: 

^ T(B:{t)^f)-r{B,{t)^f) ^ _ ^g|,cos(Amt)+^gg-sin(Amt) 
'''^''^ ^ r(5,(t) ^/) + r(5,(t) ->/) cosh(Art/2) + ^Arsinh(Art/2)' ^ ^ 

Here the mass and width difference Am and AF can be found in (|]) and A^p, A^p and have 
been defined in (P^. We have set the small parameter a in to zero and will continue to do so. The 
final states B^ DfDj, if) r]^'\ -ipfo or Xco<P in Table |T] are CP eigenstates. Their CP eigenvalue r]f 
reads Vd+d^ = = = +1 ^^'^ V^fo =Vxco^ = ^^^^ (EZP ^^^^ ^^'^ ^^'^^ 



/,^ ^ Vf sin0 sin(Amt) 

"^^^ ^ cosh(Art/2) -77/|cos0|sinh(|Ar|t/2)' ^ ^ 

Since AF and cos (p have the same sign (see (P9[)) we could replace these quantities by their absolute 



values in the denominator of (p9[). This displays that the ambiguity in the sign of cos0 cannot 
be removed by measuring acp. Its measurement determines sin0 and leaves us with a two-fold 
ambiguity in (p. Then we still do not know whether the heavier or lighter mass eigenstate is shorter- 
lived. The resolution of this ambiguity will be discussed in Section |5|. 



4.2 The CP-violating Observables of J/tl^ (f) and D*+D*- 

The situation in the decay Bs J /ip (p, which is very promising for fi-physics experiments at 
hadron machines because of its nice experimental signature, is a bit more involved than in the case 
of the pseudoscalar-pseudoscalar modes Bg DfD^ and J/tpr]^'\ since the final state is an ad- 
mixture of different CP eigenstates. In order to disentangle them, we have to make use of the angular 
distribution of the decay products of the decay chain Bg J/iJj[-^ /+/^] 0[— > K^K^], which can 
be found in [|rT],|T^. In that paper, also appropriate weighting functions are given to extract the 
observables of the angular distribution in an efficient way from the experimental data. For an ini- 
tially, i.e. at time t = 0, present i?s -meson, the time-dependent angular distribution can be written 
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generically as 



/(e,$,vi/;t) = 5]o(^)(t)^(^)(e,<i>,vi/), 



(60) 



where we have denoted the angles describing the kinematics of the decay products of J/^ ^ l^l^ 
and — i> K^K~ by 9, $ and ^. The observables 0'^^\t) describing the time evolution of the 
angular distribution ( |60| ) can be expressed in terms of real or imaginary parts of certain bilinear 
combinations of decay amplitudes. In the case of decays into two vector mesons, such as Bg 
J/ipcf), it is convenient to introduce linear polarization amplitudes Ao{t), A\\{t) and A±(t) [p]]. 
Whereas A±{t) describes a CP-odd final-state configuration, both Ao{t) and Ay (t) correspond to CP- 
even final-state configurations. The observables O*^*") (t) of the corresponding angular distribution are 
given by 

\Af{t)f with /G{0,||,±}, (61) 
as well as by the interference terms 

Re{A*{t)A\\{t)} and lm{A}{t)A^{t)} with /g{0, ||}. (62) 

For our consideration, the time evolution of these observables plays a crucial role. In the case of the 
observables (|6T|), which correspond to "ordinary" decay rates, we obtain 



\Ao{t)\' = |Ao(0)|V 



-rt 



cosh 



Art 



cos0|smh + sm0sm(Amt) 



(63) 



l^iiWP 



\Ai 



, Art 

cosh — cos ( 



smh + sm0 sm(Amt) 



(64) 



\Mt)\' = \A± 



cosh 



Art 



cos 01 smh — sm0sm(Amt) 



whereas we have in the case of the interference terms (p4): 

Re{A;{t)A^^{t)} = |Ao(0)||A||(0)| cos(52 - 5i) e"^* 



(65) 



X 



, ^rt , ,, . , lArit . , . ' 

cosh — I cos 01 smh + sm0 sm(Amt) 



(66) 



lm{Ar^{t)A^{t)} = |A||(0)||Ax(0)|e- 



-rt 



sin^i cos(Amt) — cos^i cos0 sin(Amt) — cos5i sin0 sinh 



Art 



(67) 



Im{A*(t)Ax(t)} = \Aom\MO)\e 



-rt 



sin52 cos(Amt) — cos52 cos0 sin(Amt) — cos52 sin0 sinh 



Art 



(68) 
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In 6i and S2 denote CP-conserving strong phases, which are defined as follows [lTl|,p^: 

5i = arg{A,|(0)*Ax(0)}, 62 = arg{Ao(0)*Ax(0)}. (69) 

The time evolutions generalize those given in [[TT], to the case of a sizeable Bg-Bg 

mixing phase (p to cover the pursued case of new physics. A further generalization taking into 
account also the small penguin contributions can be found in . It should be emphasized that new 
physics manifests itself only in the observables 0^''\t), while the g^''\Q, $, ^)'s are not affected. 

We may use the same angles 0, $ and to describe the kinematics of the decay products of the 
CP-conjugate transition Bg Jf-ipcp. Consequently, we have 



(70) 



Within this formalism, CP transformations relating Bg [J/ (j)]f to Bg ^ [J /^p (f)]f (/ G {0, || , ±}) 
are taken into account in the expressions for the 0^''\t) and 0^''\t), and do not affect the form of 
the 5f('^)(0, $, \1/). Therefore the same functions (/(''^(G, $, ^) are present in ( |60| ) and ( [70| ) (see 



also 



). The CP-conjugate observables O (t) take the following form 



I A 



, ^rt , ,, . , lArit . , . ' 

cosh — |cos0|smh — sm0sm(Amt) 



(71) 



i^iiwr 



■ Art |Ar|t 

cosh — cos d)\ smh 



sin0 sin(Amt) 



(72) 



1 ^rt , ,, . , lArit . , . ' 

cosh + I cos 01 smh + sm(/) sm(Amt) 



(73) 



Re{Al{t)A\\{t)] = |Ao(0)| |v4||(0)| cos(52-5i)e 



-rt 



X 



1 ^rt , ,, . , lAru . , . 

cosh — I cos 01 smh — sm0sm(Amt) 



(74) 



\m{A:{t)A^{t)} = |A||(0)||Ax(0)|e 



-rt 



X 



— sin^i cos(Amt) + cos^i cos0 sin(Amt) — cos5i sin0 sinh 



Art 



(75) 



\m{Al{t)A^{t)] = \A,m\A±me 



X 



— sin52 cos(Amt) + cos52 cos0 sin(Amt) — cos52 sin0 sinh 



Art 



(76) 
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Note that one can determine sin 5i 2, cos(5i — ^2), sin (p, cos 5i cos (p. Am and | Ar| from (|531)-([7^. 
Using cos(52 — ^i) = cos 5i cos 82 + sin 5i sin §2 in ( l66| ) and ( l74l ) one realizes that these equations 
are invariant, if the signs of cos 0, AF, and cos 5i 2 are flipped simultaneously. Hence an overall 
two-fold sign ambiguity persists and the sign of cos (p remains undetermined. 

The time evolution of the full three-angle distribution of the products of the decay chain Bs 
J/ip[-^ /+/^] 0[— > K^K^] provides many interesting CP-violating observables [[T2|, The 
expressions for three-angle angular distributions can be obtained by inserting into Eqs. (64) 



and (70) of [[T|]. 

The situation is considerably simplified in the case of the one-angle distribution, which takes the 
following form [^^: 



^ oc (|Ao(t)P + |A||(t)n I (1 + cos^ e) + \A^it)\^ I sin^ 9 . 



(77) 



Here 9 describes the angle between the decay direction of the and the z axis in the J/ip rest 
frame; the z axis is perpendicular to the decay plane of ^ K^K^ . With the help of this one-angle 
distribution, the observables |ylo(t)P + |A||(t)p and |ylx(t)p, as well as their CP conjugates, can be 
determined. They provide the following CP asymmetries: 







\Ao{t)\^ + \A\m\ 


— sine/) sin (Ami) 


'\Mt)?+\A\m\ 


+ 


|^(t)P+|A|i(t)|2 


cosh(Art/2) - 1 cos 01 sinh(|Ar|t/2) 






sin0 sin(Amt) 



(78) 



(79) 



\A^{t)\^ + \A^{t)\^ cosh(Art/2) + |cos0|sinh(|Ar|t/2)' 

In contrast to these CP-violating observables, untagged data samples are sufficient to determine the 
following quantities: 



Ao{t)\' + \A\\{t)\' + \Ao{t)\' + \A\\{t)\ 



I A 



-rt 



Art . |Ar|t 

cosn cos (h\ smn 



(80) 



|A^(t)P + |A^(t)P 



2 U 



cosn h cos smn 



(81) 



Since is tiny in the Standard Model, a striking signal of new-physics contributions to Bg-Bg mixing 
would be provided by a sizeable sin either from a fit of the tagged observables ( |63| ) - (|68|), ( |7T| ) - 
([7^, or from the CP-violating asymmetries in (^P[), ( fT^ ) and ([7^, or if the untagged observables 
( |S^ ) and ( PT| ) should depend on two exponentials. Note that in ( P^ ) the coefficient of sinh(| Ar| t/2) 
is always negative. Phrased differently, the coefficient of the exponential exp(— (r + | Ar|/2)t) with 
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the larger rate is always larger than the coefficient of exp(— (F — | Ar|/2)t). In ( |RT| ) the situation is 
reversed. This feature can be used as an experimental consistency check, once AF 7^ is established. 

Let us finally note that the formalism developed in this subsection applies also to the mode Bs 
D*^ D*^, where the subsequent decay of the D*^-mesons is predominantly electromagnetic, i.e. 
D*^ — > Df'j. The corresponding angular distribution can be found in [1TT|,[T2[]. The analysis of this 
decay requires the capability to detect photons and appears to be considerably more challenging than 
that of Bg ^ J (p, which is one of the "gold-plated" channels for 5-physics experiments at hadron 
machines. Higher Dg resonances exhibiting all-charged final states, for instance Dsi{2536)^ 
D*^[-^ Dir^] K, may be more promising in this respect [|44l]. If photon detection is not possible, 
one can still distinguish -D*^'s from Df's through the energy smearing associated with the escaped 
photon [p6|]. Then one can use the lifetime method introduced in sect. |3.3.2| to find the CP-odd 



fraction x (oc |A^(0)|2) and the CP-even fraction 1 - x (oc |Ao(0)p + \A\\{0)\'^) of the D*+D*- data 
sample through (p7[). If x 7^ 1/2 there are still non-vanishing CP asymmetries, although they are 
diluted by 1 — 2x. The corresponding formula for the CP asymmetry of this weighted average of 
CP-even and CP-odd final states can readily be obtained from (p3|)-(|^ and (J7T])-(J73|): 



r{B,{t) D*+D*-) + r(5,(t) ^ D*+D*- 



2x) sin0 sin(Amt) ^^^^ 



cosh( Art/2) - (l-2x) |cos0|sinh(|Ar|t/2)' 



The same procedure can be done with the DfD*^ data sample or any other of the decay modes in 
Table ^ 

A complete angular analysis for the three-body decays in Table [T] is more involved than the 
analysis for Bg ixf). For example in Bg — > tjjKsKs, the Ks pair does not necessarily come 
from a vector resonance and could be in an S- or D-wave or even have a larger angular momentum. 
In such cases one might restrict oneself to a one-angle transversity analysis of [ p3| ] or even satisfy 
oneself with the diluted asymmetries in (F 



5 The Unambiguous Determination of 4> 

While sin </> can be measured by conventional methods, this section shows that even sign{cos cp) can 
be determined. That determination is important for various reasons. It is not only necessary for 
a complete extraction of magnitude and phase of the new physics contributions to Bg-Bg mixing, 
(p must also be known to extract the CKM angle 7 from Bg DfK^. Even if sin0 is found to 
be consistent with zero, the determination of sign(cos(f)) is necessary to distinguish the Standard 
Model prediction cos0 ^ 1 from coscp ~ —1. In the advent of new physics, sign{cos(f)) completes 
our knowledge about (p. There are several methods to extract cos 0. 



24 



In Pursuit of New Physics with Bg Decays 



Method 1: The previous section revealed that angular correlation studies of Bg — ijjcf) determine 

cos 6i cos 0. (83) 

Once sign(cos5i) is known, sign{cos(l)) follows immediately. The former can be deduced from 
theory, once first-principle calculations of 5i have progressed sufficiently [^]. Alternatively, one 
can infer sign(cos5i) from their SU(3) counterparts occurring in B^ iIjK*[-^ 7r^Ks],'4'P^,'4'^ 
decays [denoted by sign(cos 5i)], as follows: 

The angular correlations of those Bd modes are sensitive to [OHBI] 



cos Si cos 2/?. 

By applying the SU(3) relation 

sign{cos6i) = sign{cos6i), 

the relative sign between cos 2/3 and cos can be determined, but not yet the absolute sign of cos 0. 
That absolute sign can be determined, since there are methods which extract the B^-Bd mixing 
phase 2/3 unambiguously, even in the presence of new physics [pTj-pT]]. In the absence of new 
physics, (3 equals the angle [3 of the CKM unitarity triangle. In Ref. [Q, basically the same approach 
was used to determine the sign of cos 2/3. However, in that paper it was assumed that is negligibly 
small, as in the Standard Model. On the other hand, in method 1 we assume that 2(3 is known 
unambiguously, allowing the determination of cos0. Using a theoretical input p6[ ] to determine 
sign(cos Si) as noted above, the angular distribution of the Bd — > J/tp{^ l~^l~)K*^{^ tt^K^) 



decay products considered in Ref. [ p2| ] also allows an unambiguous determination of 2/3 in the 
presence of 7^ 0. 

Method 2: Consider certain three- (or n-) body modes / that can be fed from both a Bs and 

a Bg, and where the Bs -decay amplitude is a sum over a non-resonant contribution and several 
contributions via quasi two-body modes. The strong phase variation can be modelled by Breit- 
Wigners and is known, so that cos can be extracted. Such a method was suggested in determining 
cos 2a and cos 2(3 in Bd decays [|5T|]. 

An additional method can be found elsewhere [|53|]. 



6 Conclusions 

In this paper we have addressed the experimental signatures of a non-vanishing CP- violating phase 
in the Bg-Bg mixing amplitude. Since is negligibly small in the Standard Model, but sizeable in 
many of its extensions, it provides an excellent ground for the search of new physics. We have dis- 
cussed the determination of from both untagged and tagged decays in CP-specific Bg decay modes 
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triggered by the dominant quark level decays b ccs and h cud. From lifetime measurements 
in these modes one can find the product of cos (f) and the width difference AF in the Bs system. 
The previously proposed methods to separately determine |AF| and | cos0| from untagged decay 
modes require two-exponential fits to the time evolution of either flavour- specific or CP-specific de- 
cay modes. In both cases terms of order (AF)^ must be experimentally resolved, which requires a 
substantially higher statistics than needed to measure AF cos 0. We have proposed a new method 
to measure |AF| and | cos0|, which only requires lifetime fits to the collected data samples with 
double-charm final states. This method does not require sensitivity to C((AF)^) terms. It is based 
on the observation that the measurement of AF from branching ratios discussed in [ [T6[ ] and per- 
formed in [ ^7] ] is almost unaffected by new physics. These branching ratios and AF cos0 obtained 
from the lifetime fits allow one to solve for | AF| and | cos 0| . In this context we have stressed that the 
lifetime measurements also allow one to determine the size of the CP-even and CP-odd components 
of Dl^D*~ and DfDl^ final states. This is relevant for experiments which cannot detect photons 
well enough and therefore cannot separate these components with angular analyses. We have further 
mentioned that a non-zero phase </> leads to tiny Ami oscillations in untagged data samples. This 
implies that in principle the measurement of CP violation in mixing from flavour-specific decays 
does not require tagging. 

For the tagged analyses we have generalized the formulae for the CP asymmetries to the case of 
a non-zero (p. Here we have discussed in detail the expressions needed for the angular analysis in 
Bg ilxp decays or other final states composed of two vector particles. Finally we have shown how 
the discrete ambiguities in encountered with the measurements of | cos 0| and sin can be resolved 
and can be determined unambiguously. This is important, even if sin is found to be consistent 
with zero, because it distinguishes the Standard Model case ~ from the case ~ tt. If there 
are new particles which couple to quarks with the same CKM elements as W bosons, there can be 
new contributions to the Bg-Bg mixing amplitude with larger magnitude, but opposite sign than the 
Standard Model box diagram. In this case one encounters ~ tt. This situation can occur in multi- 
Higgs doublet models and in supersymmetric models with flavour universality. From a measurement 
of Am alone the contributions from the Standard Model and from new physics to the Bg-Bg mixing 
amplitude cannot be separated. The new contribution can only be determined by combining the 
measurements of Am and 0. Consider, for example, that Am is measured in agreement with the 
Standard Model prediction: the new physics contribution to Bg-Bg mixing then varies between 
and twice the Standard Model prediction, if is varied between and ±7r. 
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